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Abstract. We discuss the universal critical behavior in (2+l)-flavor QCD by analyzing lattice data from improved staggered 
fermions generated by the HotQCD Collaboration. We present recent results from two different lattice discretizations and 
various lattice spacings (N t = 6,8, 12) at fixed physical strange quark mass (m s ) but varying light quark mass (mi). We find 
that the chiral order-parameter, i.e. the chiral condensate, shows the expected universal scaling that is associated with the 
critical point in the chiral limit already for light quark masses m//m s <0.05. From an analysis of the disconnected chiral 
susceptibility we estimate a preliminary value of the QCD transition temperature. 
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INTRODUCTION 

For almost 3 decades, lattice QCD is providing us with 
valuable information on QCD thermodynamics. Only re- 
cently, however, the numerical calculations arrive at a 
precision level, where controlled continuum extrapola- 
tions become possible. Moreover, our understanding of 
lattice cut-off effects resulting from the discretization of 
the QCD Lagrangian developed over the years. 

We will present here results from two improved stag- 
gered actions: asqtad, which has been widely used for 
finite temperature calculations in the past [1] and the 
recently introduced highly improved staggered quark 
(HISQ) action [2]. Apart from different treatments in the 
gauge part, these actions mainly differ in their strate- 
gies for the reduction of the flavor symmetry breaking, 
which is inherent to the staggered fermion formulation. 
The HISQ action has been shown to offer the most supe- 
rior degree of improvement in the flavor symmetry sec- 
tor [2, 3], In the following we will refer to the action 
that incorporates the HISQ improvement program in the 
fermion sector together with a tree-level Symanzik im- 
provement in the gauge part as HISQ/tree. 

At sufficiently high temperature, QCD undergoes a 
transition form hadronic to quark degrees of freedom, 
which is a true phase transition only in the limit of mass- 
less quarks. This transition, driven by the restauration 
of the chiral symmetry, is of second order in the chi- 
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ral limit. All thermodynamic quantities that can be de- 
rived as derivatives from the QCD partition function are 
thus subject to a universal critical behavior. In the fol- 
lowing we will analyze the universal critical behavior of 
the chiral condensate, which is a true order-parameter in 
the chiral limit. Furthermore we estimate the transition 
temperature from the disconnected part of the chiral sus- 
ceptibility. 

The lattice data has been generated by the HotQCD 
collaboration [4] with the above mentioned actions on 
lattices with temporal extent of N T = 6,8 and for the 
asqtad action with N T = 12. This translates into lattice 
spacings a of roughly a = 0.16,0.12 and 0.08 fm (at T = 
200 MeV). Throughout all calculations the strange quark 
mass m s was fixed to its physical value, whereas we have 
varied the light quark mass m; in the range mi/m s = 
0.05 — 0.2. The physical value of light to strange quark 
mass is approximately given by mi/m s = 0.037. 

UNIVERSAL SCALING 

The QCD phase transition in the chiral limit of 2-flavor 
QCD is expected to be in the 0(4) universality class [5]. 
For a sufficiently large strange quark mass this is ex- 
pected to be the same in (2+l)-flavor QCD. The stag- 
gered fermion formulation on the lattice brakes chiral 
symmetry. For each fixed lattice spacing one is only 
left with a 0(2) symmetry. The 0(4) symmetry will, of 
course, be restored in the continuum limit. 

In order to analyze the universal critical behavior that 
arises in the chiral limit we have to arrange the free en- 
ergy (and its derivatives) in terms of the relevant scaling 
variables. In general one separates two kinds of contri- 




FIGURE 1. Left: Chiral order-parameter as defined in Eq. (3) for different quark masses, piotted as a function of temperature. 
Simulations have been performed with the asqtad action obtained on N t = 8 lattices. Right: Disconnected part of the chiral 
susceptibility for ni[/m s = 0.05. Compared are results from two different lattice actions and different lattice spacings. 



butions to the free energy density, a part (f s ) that will 
generate singularities in higher order derivatives of the 
partition function and a regular part (f r ), we define 

f(T,mi,m s ) = f s (T,mi,m s )+f r (T,mi,m s ) . (1) 

Although f s depends on many parameters of the QCD 
Lagrangian, the universal behavior can be expressed in 
only two relevant scaling variables t and h, that con- 
trol deviations from criticality at (t,h) = (0,0). The scal- 
ing variables are up to some normalization factors to, ho 
given by the reduced temperature t = (T — T c )/ (toT c ) and 
the symmetry braking parameter, which we will define 
here as h = mi/(hotn s ). The universal critical behavior 
of the order parameter, M ~ df/dm/ , is controlled by a 
scaling function fc that arises from the singular part of 
the free energy density after taking a derivative with re- 
spect to the light quark mass. We make use of the fact that 
f s is a homogeneous function of its arguments, /g (and 
f s ) can thus be written as a function of the single scaling 
variable z = t/h l / s P, where 8 and j3 are critical expo- 
nents of the three dimensional 0(N) universality class. 
In addition we also consider the leading oder contribu- 
tion coming from the regular term of the free energy. We 
finally arrive at 

M = h l l s f G {z)+a,th + bih (2) 

The scaling function fc is well known from spin-model 
simulations and easily accessible via, e.g., the implicit 
parameterization given in [6]. On the right hand site 
of Eq. (2) we are left with 5 unknown parameters, 
to,ho,T c ,a t ,bi, which have to be determined through a 
fit to the lattice data. Out of these, T c is the only one 
that does not depend on the particular definition of the 
order-parameter M. All of them are, however, cut-off de- 



pendent and needed to be extrapolated to the continuum 
limit. 

On the lattice the chiral order parameter (chiral con- 
densate) is always finite but contains contributions that 
diverge in the continuum limit, i.e. it requires renormal- 
ization and in particular an additive and multiplicative 
renormalization. In order to remove at least the mul- 
tiplicative renomalization factor already on the lattice 
we consider here the following definition of the order- 
parameter 

M^m^y^jT 4 . (3) 

In Fig. 1 (left) we show the chiral oder-parameter as de- 
fined in Eq. (3) from simulations with the asqtad action 
obtained on lattices with temporal extend N T — 8 and 
light quark masses ranging from mi/m s = 0.05 — 0.2. 
Also shown are fits to the 0(2) (solid curves) and 0(A) 
(dashed curves) scaling functions. In the fits we have 
omitted the data set with the largest quark mass ;«/ jm s 
0.2 For masses to/ /m s <Q.l the fit describes the lattice 
data rather good in the entire temperature range we have 
considered. This coincides with finding obtained with the 
p4 action on N z — A [7] and N z — 8 [8]. Using the fits, we 
can extrapolate the chiral order-parameter to the physical 
mass (to//to s = 0.037) and the chiral limit, as also indi- 
cated in the Fig. 1 (left). The main difference between 
the fits manifests in the chiral limit. We obtain values 
for pranging from T c « (165 — 170) MeV (obtained for 
N T = 8 with the asqtad action). 

THE CRITICAL TEMPERATURE 

It is the ultimate goal and a valuable input for many 
phenomenological calculations, e.g. in heavy ion phe- 
nomenology, to obtain a continuum extrapolated value 




FIGURE 2. Transition temperature as obtained from the 
peak position of the disconnected part of the chiral suscepti- 
bility for different light quark masses and lattice spacings from 
simulations with the asqtad action. Also show is a combined 
chiral and continuum fit with the Ansatz given in Eq. 5. 



for T c . The chiral susceptibility is traditionally used to 
determine the QCD transition temperature. This quantity 
exhibits a peak as plotted as a function of temperature, 
which position can be associated with the QCD transi- 
tion. Note that at nonzero light quark mass the QCD tran- 
sition is a smooth crossover rather then a genuine phase 
transition of second Order. In this case the correspond- 
ing pseudo-critical temperature is not unique; it depends 
on the observable used to define it. All definitions will, 
however, lead to the same value for T c in the chiral and 
continuum limit. 

The chiral susceptibility which has a connected and 
disconnected part is defined as the mass derivative of the 
order-parameter 



T d ._ . 

X^TTj— (YY)l = Zdis + Zcon 

V ami 



(4) 



In Fig. 1 (right) we plot the disconnected part of the 
chiral susceptibility as obtianed for light quark mass 
mi/m s = 0.05 but different lattice spacings. We com- 
pare results from two different lattice actions: asqtad and 
HISQ/tree. We find that the cut-off effects are visible in 
the peak hight and peak position. They are considerably 
smaller in case of the HISQ/tree action. 

From the scaling function for the chiral condensate, 
Eq. (2), the scaling of the chiral susceptibility can be 
easily extracted. It is governed by the function f G (z) = 
dfa I dz. To obtain a consistent picture it is desirable to 
perform a similar analysis as presented for the order- 
parameter in the previous section. This should lead to the 
same non-universal normalizations of the scaling fields 
as well as to the same value for T c in the chiral limit. At 
the moment, however, we restrict ourselves to an analysis 
of the peak positions T p of the disconnected part x&s- 



From Eq. (2) we obtain the fit Ansatz 

T p (m h N % ) = T c + b(m l /m s ) d + cN. 
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(5) 



Here we allow for a generic cut-off dependence of the 
form 1/N% in addition to the mass dependence coming 
from the scaling Ansatz. We thus have three fit parame- 
ters, T c , b, c. The exponent d is fixed by universal scaling 
asd = l/j85. 

In Fig. 2 we show the peak positions T p from 
as obtained by simulations with the asqtad action for 
different quark masses and lattice spacings. The values 
and their errors have been estimated by two different 
asymmetric fits to the peak [9]. In addition we show a 
combined chiral and continuum extrapolation based on 
Eq. (5). The fit works reasonably well and we obtain 
as a preliminary result for the transition temperature at 
physical masses: 



T p (0.037m,, oo) = 164(6) MeV. 



(6) 



Here the error summarizes our estimate of the statistical 
and systematical error. 
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